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ABSTRACT 


A 

Metal-ultrasmall  superlattice  junctions  have  been 
analysed  by  using  a  consistent  tight-binding  method.  Hetero¬ 
junctions  and  metal-semiconductor  interfaces  have  been  discussed 
independently  by  assuming  that  each  semiconductor  superlayer 
has  more  than  six  layers. 

We  present  results  for  the  electronic  structure  of 
( 100 ) -semiconductor  hetero junctions ;  in  particular,  we  consider 
the  (lOO)-GaAs-AlAs,  GaSb-InAs  and  InAs-GaSb  interfaces.  For 
GaAs-AlAs  we  have  found  that  the  valence  band  offset  is  0.41  eV, 
and  that  there  is  an  interface  dipole  of  0.13  eV  (for  no  inter¬ 
face  dipole,  the  valence  band  offset  would  be  of  0.28  eV).  Our 
results  also  show  that  the  two  GaSb-InAs  and  InAs-GaSb  interfaces 
present  different  band-offset  discontinuities. 

We  have  also  analysed  the  initial  stages  of  the 
Schottky-barrier  formation  for  A1  on  an  As-terminated  GaAs-(lOO) 
face.  Two  cases  have  been  analysed:  (i)  firstly,  we  have 
considered  a  half  monolayer  of  Al,  with  these  atoms  located  in 
the  ideal  position  corresponding  to  a  Ga-atom  of  the  ideal 
continued  GaAs  crystal;  (ii)  in  the  second  case,  we  have 
considered  a  full  monolayer  with  the  Al  atoms  forming  a  (100)- 
face  of  an  Al  crystal. iOur  results  show  that  the  Schottky  barrier 
is  formed  for  a  monol^ybr,  and  that  the  barrier  height  is  0.55 
eV,  similar  to  the  value  rpund  for  a  metal-GaAs ( 110 )  junction. 

The  consistency  found  for  the  GaAs-AlAs ( 100 )  and  the 
Al-GaAs(lOO)  interfaces  defines  completely  the  barrier  heights 
associated  with  the  Al-GaAs-AlAs ( 100 )  junction. 
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1.  INTRODUCTION 


The  mechanism  of  the  Schottky  barrier  formation  is  a 
matter  of  current  controversy.  Different  models  have  been 
proposed,  and  it  is  clear  that  none  has  a  universal  acceptance 
since  no  single  model  can  explain  all  the  experimental  evidence. 
This  is  related  to  the  different  kind  of  interfaces  appearing 
in  a  metal-semiconductor  junction;  according  with  the 
experimental  evidence  there  are:  (i)  abrupt  interfaces  with 
a  well  defined  separation  between  the  metal  and  the  semiconductor 
and  (ii)  reactive  interfaces  3'6  with  the  metal  and  the  semi¬ 
conductor  interdif fusing  and/or  forming  a  new  chemical 
compound  at  the  interface. 

Two  different  models,  related  to  the  two  kind  of  inter¬ 
faces  commented  on  above,  appear  to  attract  fairly  wide  support 
in  the  scientific  community:  (i)  The  intrinsic  interface  states 
model  7-11 ,  which  attributes  the  Schottky  barrier  height  to  the 
Fermi  level  pinning  by  the  metal  induced  density  of  states  near 
the  semiconductor  mid-gap.  (ii)  The  defect  model  12-16 f  which 
assumes  that  the  Schottky  barrier  formation  results  from  the 
creation  of  a  significant  number  of  defects  at  or  near  the 
interface.  Obviously,  the  defect  model  is  assumed  to  be  valid 
for  reactive  interfaces,  while  the  intrinsic  interface  states 
model  seems  to  be  related  to  abrupt  interfaces. 

In  this  work,  we  have  started  an  investigation  of  the 
Schottky  barrier  formation  in  metal-ultrasmall  superlattices , 
for  cases  for  which  there  is  an  abrupt  interface  between  the 
metal  and  the  superlattice . 

Different  theoretical  approaches  have  been  used  to 
analyse  those  ideal  abrupt  interfaces.  It  is  worth  mentioning 
three  different  approaches:  (i)  in  the  simplest  case,  the 
mechanism  of  the  Schottky  barrier  formation  8,11  is  based  on  the 
concept  of  the  charge  neutrality  level.  Here,  when  a  semiconduc¬ 
tor  and  a  metal  are  intimately  coupled,  the  semiconductor 
density  of  surface  states  broadens  but  keeps  its  centre  of 
gravity  almost  constant.  Hence,  the  interface's  Fermi  level 
almost  coincides  with  the  Fermi  level  for  the  free  surface.  Then, 
the  barrier  height  is  determined  by  equalising  both  the  metal 


Fermi  level  and  the  semiconductor  charge  neutrality  level. 

(ii)  In  the  most  elaborate  kind  of  calculation  ^ ,  the  barrier 
height  is  determined  by  means  of  a  self-consistent  calculation 
using  pseudopotential  theory  and  the  local  density  formalism. 

(iii)  In  an  intermediate  approach  9»10f  we  have  used  recently 
a  consistent  calculation  for  a  tight-binding  model.  The  idea 
behind  this  approach  is  to  extend  the  charge  neutrality  level 
concept  8  to  a  more  sophisticated  calculation,  which  takes  into 
account  the  small  variations  that  interface  dipoles  introduce 
in  the  "rule  of  thumb"  associated  with  the  charge  neutrality 
approach. 

In  the  work  presented  in  this  project  we  have  followed 
this  last  method:  the  "consistent  tight-binding  approach",  a 
method  that  has  been  applied  successfully  to  metal-semiconductor 
junctions  9,10.  on  the  other  hand,  we  have  analysed  the  metal- 
ultrasmall  superlattice  junctions  in  several  steps:  (i)  in  a 
first  step,  we  have  discussed  the  properties  of  free-surface 
semiconductors.  The  purpose  of  this  work  is  to  analyse  the 
consistency  of  the  tight-binding  method  in  a  simple  case,  and 
to  find  out  how  the  surface  perturbation  penetrates  into  the 
crystal,  (ii)  In  a  second  step,  we  have  considered  the  semicon¬ 
ductor-semiconductor  interface,  assuming  it  to  be  decoupled  from 
the  metal-semiconductor  junction.  The  idea  behind  this  approach 
is  that  the  metal-superlattice  junction  properties  are  determined 
by  the  independent  semiconductor-semiconductor  and  metal-semi¬ 
conductor  interfaces  (we  only  consider  in  this  work  a 
configuration  with  the  metal-semiconductor  interface  parallel 
to  the  semiconductor-semiconductor  interface  of  the  superlattice) 
As  our  results  show,  this  is  a  valid  assumption  if  each  semi¬ 
conductor  superlayer  has  more  than  six  semiconductor  layers. 

(iii)  In  a  final  step,  we  consider  the  metal-semiconductor 
junction  and  analyse  the  Schottky  barrier  height  for  the  initial 
stages  of  the  barrier  formation. 

In  the  following  paragraphs  we  present  a  brief  discussion 
of  the  results  obtained  in  the  different  papers  collected  in 
the  last  part  of  this  report.  In  §2,  paper  No.l:  "Anion-induced 
surface  states  for  the  ideal  (100) -faces  of  GaAs ,  AlAs  and 
GaSb",  by  G.Platero,  J.S£nchez-Dehesa,  C.Tejedor,  F. Flores  and 


A.Mufioz,  is  presented.  In  S3,  we  discuss  papers  No. 2  and  3: 

"Electronic  structure  of  ( 100) -semiconductor  hetero¬ 
junctions",  by  G.Platero,  J.S4nchez-Dehesa,  C.Tejedor  and  P. 
Flores ,  and 

"Band-offset  calculations  of  (100) -semiconductor  hetero 
junctions:  A  simple  consistent  method",  by  G.Platero,  J.Sinchez- 
Dehesa,  C.Tejedor  and  F. Flores. 

Finally,  in  §4,  we  discuss  papers  No. 4  and  5: 

"Initial  stages  of  the  Schottky  barrier  formation  for 
abrupt  covalent  interfaces",  by  G.Platero,  J. A. Verges  and  F. 
Flores,  and 

"Initial  stages  of  the  Schottky  barrier  formation  for 
A1  on  GaAs(lOO)",  by  G.Platero,  J.Durdn  and  F. Flores. 

In  §5,  we  present  our  final  comments. 

2.  PAPER  No . 1 :  "Anion-induced  surface  states  for  the  ideal 
(100) -faces  of  GaAs,  AlAs  and  GaSb" . 

In  this  paper  we  apply  our  consistent  tight-binding 
method  to  analyse  the  anion-induced  surface  states  of  the  ideal 
(100) -faces  of  GaAs,  AlAs  and  GaSb. 

Our  results  show  that  consistency  is  crucial  in  order 
to  obtain  the  surface  state  levels  and  the  Fermi  energy.  For 
example,  for  GaAs  we  have  found  that  the  Fermi  level  is  located 
at  0.33  eV  above  the  valence  band  top;  had  we  not  applied  our 
consistent  procedure,  we  would  have  obtained  the  Fermi  level  at 
0.77  eV  above  the  valence  band  top. 

On  the  other  hand,  the  results  of  paper  No.l  show  that 
the  surface  perturbation  penetrates  three  layers  in  the  semi¬ 
conductor,  in  such  a  way  that  the  fourth  layer  is  practically 
unaltered  by  the  surface.  This  is  an  important  result  that  will 
be  used  later  on,  in  order  to  analyse  semiconductor  interfaces. 

3.  PAPERS  No. 2  AND  3:  "Electronic  structure  of  (100 ) -semiconduc 
tor  hetero junctions"  and  "Band-offset  calculations:  A  simple 
consistent  method". 

In  these  papers  we  apply  the  consistent  tight-binding 
method  to  the  calculation  of  the  band-offsets  for  the  (100)- 
GaAs-AlAs,  GaSb-InAs  and  InAs-GaSb  hetero junctions .  In  our 
method  we  have  included  spin-orbit  coupling  in  the  parameters 


of  the  different  semiconductor's  Hamiltonians  and  used  sp  s* 
hybrids  in  each  atom. 

Each  hetero junction  is  analysed  by  projecting  the  whole 
structure  of  both  semiconductors  onto  the  six  layers  located 
around  the  interface,  three  layers  corresponding  to  one  semi¬ 
conductor.  This  means  that  the  whole  Hamiltonian  for  both 
crystals  is  reduced  to  an  effective  matrix-Hamiltonian  having 
60x60  elements. 

Consistency  is  introduced  by  means  of  diagonal  perturba 
taions  at  the  interface  layers.  Using  the  effective  matrix- 
Hamiltonian  we  calculate  the  electronic  charge  at  each  layer, 
and  define  consistency  relating  the  diagonal  perturbations  with 
the  induced  charge  through  a  kind  of  electrostatic  (or  Hartree) 
equation. 

For  GaAs-AlAs,  we  have  found  that  the  valence  band 
offset  is  0.41  eV,  and  that  there  is  an  interface  dipole  of 
0.13  eV  (for  no  interface  dipole,  the  valence  band  offset  would 
be  of  0.28  eV) . 
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Recent  experimental  information  yields  a  value  of 
0.65  eV  for  the  valence  band  offset.  It  must  be  noticed  that 
this  value  has  been  accepted  long  since  to  be  0.25.  We  think 
that  our  results  (including  an  uncertainty  of  0.1  eV  associated 
with  the  approximations  of  our  theoretical  model)  are  quite 
satisfactory. 

4.  PAPERS  No. 4  AND  5:  "Initial  stages  of  the  Schottky  barrier 
formation  for  abrupt  covalent  interfaces"  and  "Initial  stages 
of  the  Schottky  barrier  formation  for  A1  on  GaAs(100)" . 

In  these  papers  we  have  analysed  the  initial  stages  of 
the  Schottky  barrier  formation  for  A1  on  Si  (111)  and  GaAs(100). 

Si (111)  has  been  chosen,  in  a  first  step,  as  a  well- 
known  case  for  analysinq  the  initial  staqes  of  the  Schottkv 
barrier  formation.  Aqain,  we  have  applied  the  consistent  tiqht- 
bindinq  method  and  analysed  the  evolution  of  the  electronic 
density  of  surface  states  in  the  semiconductor  qap  as  a  function 
of  the  deposition  of  Al  on  Si.  Two  cases  have  been  considered: 

(i)  a  third  of  monolayer,  and  (ii)  a  full  monolayer.  It  is 


well-known  that  for  a  third  of  monolayer  there  is  an  energy 
gap  in  the  density  of  states,  in  such  a  way  that  the  barrier 
height  does  not  appear  as  completely  formed.  Our  analysis  for 
a  full  monolayer  have  shown  that:  (i)  the  barrier  is  completely 
formed  in  this  case,  with  the  electronic  density  of  states 
filling  the  semiconductor  energy  gap,  and  (ii)  the  charge 
neutrality  level  (and  the  barrier  height)  depends  on  the 
geometrical  position  of  the  metal  atoms  on  the  silicon  surface. 

Our  analysis,  and  the  comparison  with  the  experimental  evidence, 
shows  that  the  metal  atoms,  for  a  monolayer,  are  located  in  the 
top  position. 

Finally,  we  have  considered  the  initial  stages  of  the 
Schottky  barrier  formation  for  the  case  of  A1  deposited  on  an 
anion-terminated  GaAs (100 ) -surface .  Here,  we  have  analysed  a 
half  and  a  complete  monolayer  of  Al;  for  a  half  monolayer,  an 
Al  atom  is  located  in  the  ideal  configuration  corresponding  to 
a  Ga  atom  of  the  ideal  crystal  growing  towards  the  metal.  The 
main  result  of  our  analysis  is  that,  for  a  half  monolayer,  we 
find  a  cation-like  surface  state  band  similar  to  the  one  appearing 
in  the  ideal  (100) -cation  surface.  For  a  full  monolayer,  we 
find  that  the  barrier  height  is  completely  formed,  and  that  the 
Fermi  level  is  shifted  with  respect  to  the  half  monolayer  case 
by  0.31  eV  towards  the  conduction  band.  We  conclude  that  the 
barrier  height  for  an  Al-GaAs(lOO)  face  is  0.55+0.05  eV,  similar 
to  the  value  found  for  a  metal-GaAs (11 0 )  face  (0.65±0.05  eV) . 

5.  CONCLUSIONS 

Metal-ultrasmall  superlattice  junctions  have  been 
analysed  using  a  consistent  tight-binding  method.  Energy  levels 
in  the  junction  and  the  superlattice  are  determined  by  analysing 
each  interface  which  appears  to  be  decoupled  of  each  other  if 
the  size  of  each  superlayer  is  greater  than  six  layers  of  the 
corresponding  semiconductor.  Accordingly,  hetero junctions  and 
metal-semiconductor  interfaces  have  been  discussed  independently. 
For  the  case  of  Al  on  GaAs-AlAs (100 ) ,  we  have  found  the  following 
main  results:  (i)  For  a  GaAs-AlAs  (100 )  heterojunction ,  the 
valence  band  offset  is  0.41  eV;  (ii)  for  an  Al-GaAs(lOO)  junction, 
with  Al  growing  on  an  As-terminated  crystal,  the  barrier  height 


appears  to  be  0.55  eV. 

The&e  two  vatue.6  deitne  compZeteZy  the.  AZ-GaA6-AZA6  [1 00) 
junction  and  aZZowi  tu  to  obtatn  ^  ion.  any  gtven  GaA&-AZA& 
bapen.Zatti.ee.  Thib  can  be  eabiZy  achieved  by  n.eienning  the 
eZectnonic  ZeveZ b  oi  the  bupenZattice  to  the  vaZence  band 
btnuctune  oi  the  two  ideaZ  iniinite  bemiconductoAb . 
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ANION-INDUCED  SURFACE  STATES  FOR  THE  IDEAL 
(100) -FACES  OF  GaAs,  AlAs  and  GaSb 

G.Platero,  J.S^nchez-Dehesa,  C.Tejedor,  F. Flores  and  A.Mufioz 
Departamento  de  Fisica  del  Estado  Sdlido, 
Universidad  Autdnoma,  Cantoblanco,  28049  Madrid,  Spain 


We  present  a  simple  consistent  tight-binding  method  to 
calculate  the  electronic  properties  of  ionic  semiconductor 
surfaces.  The  method  is  applied  to  calculate  the  anion-induced 
surface  states  of  the  ideal  (lOO)-faces  of  GaAs,  AlAs  and  GaSb. 

It  is  shown  that  consistency  is  crucial  to  obtain  the  surface 
state  levels  and  the  Fermi  energy,  and  that  the  surface  perturbation 
penetrates  three  layers  in  the  crystal. 
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I.  INTRODUCTION 


Many  different  theoretical  and  experimental  work  have 
been  addressed  to  the  understanding  of  the  surface  states  of 
clean  semiconductors  [l,2] .  Most  of  the  work  have  been  done  in 
covalent  surfaces  and  GaAs  [3]  .  Specific  work  on  the  (lOO)-surfaces' 
of  ionic  crystals  have  been  stimulated  by  Molecular  Beam  Epitaxy 
and  the  growing  of  superlattices. 

The  main  experimental  data  for  the  (100) -surface  has 
been  collected  in  GaAs  [4,5],  but  interesting  information  for 
other  semiconductors,  GaSb  and  InAs,  is  also  available  [6,7]. 

From  the  theoretical  point  of  view,  let  us  mention  the 
early  work  of  Appelbaum  et  al.  [8] ,  where  a  self-consistent 
solution  for  the  cation-like  ( 100 ) -surface  of  GaAs  is  presented. 
Ivanov  et  al.  I3]  have  calculated  the  cation-like  and  the  anion¬ 
like  surface  states  of  the  GaAs  (lOO)-face  by  means  of  a  tight- 
binding  approach,  but  no  kind  of  consistency  is  provided  in  this 
calculation. 

In  this  paper  we  have  addressed  the  problem  of 
calculating  the  anion-like  surface  states  of  GaAs,  AlAs  and 
GaSb.  It  is  well  known  that  the  (lOO)-faces  of  those  crystals 
present  different  reconstructions  depending  on  the  kind  of 
preparation  [4-7] .  It  has  been  suggested  [7]  that  the  different 
GaAs  surface  reconstructions  are  associated  with  a  kind  of 
metal-insulator  transition:  this  means  that  dangling-bonds  are 
paired  to  produce  bonds  having  no  density  of  states  at  the  Fermi 
level.  For  GaSb(lOO)  surfaces,  however,  reconstructions  seem  to 
be  related  with  a  three-fold  or  a  five-fold  periodicity  [7] : 
unless  the  presence  of  vacancies  is  involved,  this  means  that 
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the  electronic  properties  of  these  surfaces  must  differ 
markedly  of  the  GaAs(100)-faces. 

In  spite  of  these  comments,  we  have  only  considered 
in  this  paper  the  ideals  anion-like  (lOO)-surfaces  of  GaAs,  AlAs 
and  GaSb  [9] ;  we  consider  these  surfaces  to  be  of  interest  by 
themselves,  mainly  considering  that  ideal  crystals  (with  no 
reconstruction)  appear  in  a  superlattice  built  up  with  these 
semiconductors.  The  purpose  of  this  paper  is  to  present,  for 
these  surfaces,  a  simple  method  which  includes  the  basic 
consistent  properties  associated  with  charge  neutrality 
conditions.  By  means  of  the  self-consistent  conditions  imposed 
at  the  surface  we  also  analyse  how  deep  in  the  crystal  the 
surface  perturbation  penetrates,  and  the  Fermi  level  position, 
a  quantity  very  sensitive  to  the  charge  neutrality  conditions 
and  intimately  related  to  the  charge  neutrality  level  of  the 
semiconductor  [lO]  . 

The  calculation  presented  in  this  paper  for  the  anion¬ 
like  surfaces  has  been  done  by  using  a  simple  tight-binding 
approach  (similar  to  Ivanov  et  aL’s  calculation  [3])  ,  bat  supplemented 
by  a  kind  of  self-consistency  appropriate  for  this  model.  This 
consistency,  and  the  charge  neutrality  conditions,  has  been 
proved  to  be  crucial  for  the  understanding  of  other  interfaces 
[ll,12].  The  main  idea  is  to  introduce  diagonal  perturbations 
at  several  layers  (the  layers  affected  are  determined  specifically 
at  each  case)  at  the  interface.  These  perturbations  are 
determined  by  introducing  self-consistency  between  the  potentials 
and  the  surface  charges. 

In  §2  we  discuss  our  model,  our  method  of  calculation 
and  how  self-consistency  is  introduced  in  the  problem.  In  S3  we 
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present  our  results  and  our  conclusions. 

II.  THE  MODEL,  METHOD  OF  CALCULATION  AND  CONSISTENCY 

We  describe  the  semiconductors  by  means  of  a  tight- 
binding  model,  using  sp^s*  hybrids  [l3]  in  each  atom  and  include 
interaction  parameters  extending  up  to  the  first  neighbours. 
Table  1  gives  the  parameters  used  in  our  calculation  for  GaAs, 
AlAs  and  GaSb.  We  assume  that,  except  for  the  diagonal 
perturbations  commented  on  above ,  at  the  interface  we  have  the 
same  parameters  as  in  the  bulk.  In  our  calculation  we  have  also 
included  the  spin-orbit  interaction  [l4]  that  appears  to  be 
important  for  Sb.  The  parameters  given  in  Table  I  are  similar 
to  the  ones  proposed  by  Vogl  et  al.[l3]  except  by  small 
corrections  introduced  to  allow  for  the  spin-orbit  interaction 
(the  parameters  given  in  Table  1  have  been  chosen  to  give  the 
right  levels  at  the  symmetric  points  T  and  X). 

The  Hamiltonian  defining  our  surface  model  is, 
therefore,  well  known  except  for  the  diagonal  perturbations 
applied  at  some  layers  of  the  interface  (see  below). 

As  regards  our  method  of  calculation,  we  follow  the 
procedure  explained  in  ref. [l5] .  In  summary,  our  method  consists 
in  looking  for  the  surface  components  of  the  Green  function  of 
the  system.  By  means  of  a  decimation  technique,  we  can  project 
the  bulk  components  of  the  Green  function  onto  a  few  layers  of 
the  interface.  In  our  actual  case,  we  have  projected  all  the 
electronic  bulk  structure  onto  the  last  three  semiconductor 
surface  layers.  This  means  that  the  whole  Hamiltonian  is  reduced 
to  an  effective  matrix  Hamiltonian  having  30x30  elements 


(including  spin  and  five  orbitals  per  layer).  This  matrix 
defines  the  effective  interactions  for  the  surface  atoms,  in 
such  a  way  that  now  we  analyse  the  surface  by  introducing  in 
this  effective  Hamiltonian  the  diagonal  perturbations  to  be 
discussed  just  below. 

Turning  our  attention  to  self-consistency,  let  us  start 
our  discussion  by  considering  the  electrostatic  potential 
associated  with  the  ideal  charges  existing  in  the  bulk  when 
assumed  to  be  extended  until  the  same  surface.  Figure  1  shows 
this  situation  after  averaging  in  the  plane  parallel  to  the 
surface.  Notice  that  for  a  cation  we  assume  to  have  (3+a)  units 
of  charge  per  atom,  and  (5-a)  units  of  charge  per  atom  for  an 
anion  excep-t  for  the  last  layer.  It  is  well  known  [16]  that  for 
a  (lOO)-ionic  face,  the  charge  of  the  last  anion  must  be  (5-y), 
this  value  being  determined  by  the  condition  of  having  a  finite 
potential  inside  and  outside  the  crystal.  Let  us  also  comment 
that  this  surface  charge  is  also  associated  with  the  surface 
state  occupancy;  for  the  (100) -anion  surface  we  find  two  surface 
states  (see  below)  with  a  total  occupancy  of  5/4  (one  is  fully 
occupied  and  the  other  one  only  in  1/4). 

However,  the  ideal  charges  shown  in  Figure  1  are  not 
consistent  with  the  potentials  of  an  ideal  semiconductor 
extending  unperturbed  up  to  the  surface.  Instead,  charges  are 
transferred,  in  our  case,  among  the  last  three  layers,  say,  1, 

2  and  3.  Now,  we  introduce  two  diagonal  perturbations,  and  V^, 
at  layers  1  and  2,  looking  for  a  kind  of  Hartree  self-consistency 
between  these  potentials  and  the  charges  transferred  among 
layers.  Thus,  calling  6n1  ,  6n9  and  <5n_  (6n1  +  6n9  +  6n«  =  0),  the 


changes  in  the -layers  charges  measured  with  respect  to  the 
ideal  values  shown  in  Figure  1  ,  we  write  down  the  following 
equations  of  consistency: 


V1  =  X  d(6n1+6n2) 


V2  =  d  &n±  .  (2) 

where  A  is  the  surface  area  per  atom. 

(Notice  that  the  charge  neutrality  condition  6n1+6n2+6n3  =  0  is 
automatically  satisfied,  up  to  a  given  accuracy,  by  imposing  the 
appropriate  surface  states  occupancies). 

Eqs.(l)  and  (2)  yield  the  consistent  equations  to  be 
used  in  order  to  determine  and  V2>  Notice  that,  on .the  other 
hand,  V*  and  V2  are  the  diagonal  perturbations  to  be  introduced 
in  the  effective  Hamiltonian  (a  30x30  matrix)  discussed  above. 

It  is  of  interest  to  note  that  a  good  approximation  to 
the  self-consistent  equations  (1)  and  (2)  can  be  obtained  by 
looking  for  values  of  and  V2 ,  say,  and  V2^\  such  that 

6n^  %  6n2  =  0.  The  idea  behind  this  approximation  is  that  a  very 
small  change  in  and  <5n2  introduces  very  strong  variations 
in  and  V2  through  the  induced  potentials  given  by  eqs.(l)  and 


III.  RESULTS  AND  COMMENTS 

Figure  2a  shows  the  local  density  of  states  for  the 
last  As-layer  of  GaAs,  as  calculated  by  using  no  consistency 
(Va  =  V2  =  0).  Figures  2b  and  2c  show  the  local  density  of  states 
projected  onto  the  p  +p„  (bridge  bond)  and  the  s+p„  (dangling  bond) 


orbitals  of  As.  From  these  figures  we  see  that  the  1/4-occupied 
surface  state  band  is  a  bridge-like  state,  while  the  occupied 
surface  state  band  is  a  dangling-bond  state,  in  agreement  with 
Ivanov  et  al's  results  [3].  Notice  that  the  Fermi  level  is 
inside  the  1/4-occupied  band  and  is  located,  for  this  non- 
consistent  case,  at  0.77  eV  above  the  valence  band  top.  It  is 
of  interest  to  comment  that  the  1/4-occupied  surface  band  affords 
a  check  to  the  accuracy  of  our  calculations,  and  to  the  number 
of  layers  to  be  used  in  order  to  obtain  a  good  consistent  result. 
In  our  actual  case,  layers  1,  2  and  3  yield  0.74  electrons 
(instead  of  0.75)  for  the  empty  number  of  states  in  the  partially 
occupied  surface  band;  had  we  used  only  two  layers,  instead  of 
three,  we  would  have  obtained  0.67  electrons.  This  shows  that 
using  three  layers  give  an  accuracy  better  than  0.01  electrons 
for  the  total  number  of  electrons  in  the  surface  band.  On  the 
other  hand,  this  discussion  shows  that  the  surface  perturbation 
penetrates  three  layers  in  the  crystal. 

Figure  3a  shows  the  local  density  of  states  for  the  As- 
layer  of  GaAs,  as  calculated  in  the  consistent  procedure.  For 
this  case,  we  have  obtained  the  following  values  for  the 
diagonal  perturbations,  and  : 

(GaAs )  =  -0.52  eV 
V ^ (GaAs)  =  0.17  eV. 

Figures  3b  and  3c  show  this  density  projected  onto  the 
bond 

bridge^ and  the  dangling-bond  orbitals,  respectively.  Note  that 
in  this  consistent  case,  Ep  =  0.33  eV,  differing  from  the 


Let  us  finally  comment  that  the  partially  occupied 
bands  for  the  three  semiconductors,  GaAs,  AlAs  and  GaSb,  are 
located  around  the  mid-gap,  and  that  their  band  widths  are 
similar,  'v  0.5  eV.  These  states  correspond  to  bridge  orbitals; 
in  this  case,  other  people  have  found  broader  bands  [3,8].  This 
is  related  to  the  kind  of  approximation  used  in  our  initial 
Hamiltonian;  we  have  taken  a  sp^s*  basis  with  interaction  up  to 
first  nearest  neighbours  only.  When  these  interactions  are 
extended  up  to  second  neighbours,  the  bridge-like  surface  band 
broadens  [3].  However,  we  can  expect,  from  the  results  for  other 
interfaces  [10,11,12"],  the  surface  Fermi  energy  and  the  charge 
neutrality  level  to  be  practically  independent  of  the  Hamiltonian 
parameters,  depending  mainly  on  the  semiconductor  electronic 
structure  for  the  bulk.  This  seems  to  be  the  case  even  for 
reconstructed  surfaces;  indeed,  the  experimental  evidence  shows 
that  the  surface  Fermi  level  for  GaAs  is  practically  independent 
of  the  actual  reconstruction:  Ep  =  0.4 ±  0.1  eV  M,  or  0 . 5  5  ±  0.1 
eV  [5]  above  the  valence  band  top.  These  values  are,  strikingly, 
in  very  good  agreement  with  the  value  obtained  above,  Ep  =  0.^3eV. 

In  conclusion,  we  have  presented  a  simple  consistent 
tight-binding  calculation  for  the  (100)-faces  of  ionic  semi¬ 
conductors.  The  method  has  been  applied  to  the  anion  faces  of 
GaAs,  AlAs  and  GaSb,  our  results  showing  that  the  surface 
perturbation  generates  three  layers  in  the  crystal.  Self-consis¬ 
tency  is  shown  to  be  crucial  in  order  to  obtain  appropriate 
values  for  the  Fermi  energy  (and  the  charge  neutrality  level),  our 
results  showing  differences  of  0.4  eV  between  the  results  for 
the  Fermi  energy  of  consistent  and  nonconsistent  calculations. 
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Figure  Cpations 

Fig.l.-  Schematic  representation  of  charges  and  electrostatic 
potentials  at  the  surface. 

Fig. 2.-  Local  density  of  states  (in  a.u.)  for  the  last  As-layer 

of  the  anion-terminated  (100)  surface  of  GaAs  with  no 

shift  of  the  surface  atomic  potentials  with  respect  to 

those  of  the  bulk,  a)  Local  density  projected  on  the  five 

atomic  orbitals  of  As.  b)  Local  density  projected  on  the 

bridge-bond  Px+Py.  c)  Local  density  projected  on  the 

dangling-bond  s+p  . 

z 

Fig. 3.-  Local  density  of  states  (in  a.u.)  for  the  last  As-layer 

of  the  anion-terminated  (100)  surface  of  GaAs  with 

consistent  potentials  V^--0.52  eV  and  V2=0.17  eV.  a)  Local 

density  projected  on  the  five  atomic  orbitals  of  As.  b) 

Local  density  projected  on  the  bridge-bond  px+py.  c)  Local 

density  projected  on  the  dangling-bond  s+p  . 
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Fig. 4.-  Local  density  of  states  (in  a.u.)  for  the  last  As-layer 

of  the  anion-terminated  (100)  surface  of  AlAs  with  consis¬ 
tent  potentials  V^=-0.17  eV  and  V2=0.21  eV.  a)  Local 
density  projected  on  the  five  atomic  orbitals  of  As.  b) 
Local  density  projected  on  the  bridge-bond  Px+Py*  c)  Local 

density  projected  on  the  dangling-bond  s+p  . 

z 

Fig. 5.-  Local  density  of  states  (in  a.u.)  for  the  last  Sb-layer 
of  the  anion-terminated  (100)  surface  of  GaSb  with 
consistent  potentials  V^-0.23  eV  and  V2=0.25  eV.  a)  Local 
density  projected  on  the  five  atomic  orbitals  of  As.  b) 
Local  density  projected  on  the  bridge-bond  p^+p^.  c)  Local 
density  projected  on  the  dangling-bond  s+p  . 


Surf,  atom 


Surf,  atom 


Rx+  Py 


ELECTRONIC  STRUCTURE  OF  (100)  SEMICONDUCTOR 
HETEROJUNCTIONS 

G.  PLATERO.  J.  SANCHEZ-DEHESA,  C.  TEJEDOR  and  F.  FLORES 

Depanamemo  de  Fisica  del  Estado  Sdlido,  Univtnidad  Auidnoma,  Camoblanco,  29049  Madrid, 
Spain 

Received  10  June  1985;  accepted  for  publication  14  June  1985 


The  electronic  structure  of  (100)  abrupt  semiconductor  heterojunctions  has  been  analysed 
by  means  of  a  self-consistent  tight-binding  method.  Band  offsets  for  the  junctions  GaAs-AIAs, 
InAs-GaSb  and  GaSb-InAs  have  been  calculated  and  compared  with  the  experimental  evidence. 


1.  Introduction 

The  extensive  use  of  synchrotron-radiation  photoemission  has  produced 
an  enormous  amount  of  experimental  information  on  semiconductor  hetero- 
junction  band  discontinuities  [1, 2].  In  spite  of  the  existence  of  a  few  self-con¬ 
sistent  calculations  to  determine  these  magnitudes  for  some  particular  cases 
[3-5],  the  theoretical  effort  is  more  concentrated  on  simpler  models  [6-8] 
suitable  for  a  large  set  of  heterojunctions.  In  general  these  models  try  to 
obtain  band  discontinuities  from  bulk  properties  including  in  some  cases  [7] 
corrections  due  to  the  charge  redistribution  at  the  interface  (IF).  The  physical 
ideas  introduced  long  ago  for  metal-semiconductor  junctions  [9]  and  semi¬ 
conductor  heterojunctions  [8]  are  properly  developed  in  this  paper  by  perfor¬ 
ming  detailed  calculations  for  heterojunctions  within  the  same  scheme  recent¬ 
ly  used  for  metal-semiconductor  junctions  [10].  The  main  contribution  to 
band  offsets  comes  from  the  alignment  of  the  charge  neutrality  points  of  both 
semiconductors,  i.e.  the  energies  up  to  which  the  density  of  states  must  be 
occupied  to  obtain  local  charge  neutrality  in  the  bulk  of  each  semiconductor. 
Afterwards,  a  non-negligible  correction  must  be  introduced  due  to  the  charge 
redistribution  in  a  few  layers  around  the  IF.  We  apply  the  method  to  GaAs- 
AIAs  and  InAs-GaSb  (100)  heterojunctions  where  some  controversy  exists 
on  the  electronic  band  structure  [11,  12]. 
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2.  Method 


We  describe  both  semiconductors  by  means  of  a  tight-binding  method, 
including  only  first-nearest-neighbour  interaction.  Self-consistency  is  intro¬ 
duced  in  the  calculation  by  means  of  different  diagonal  perturbations,  Vh  at 
five  layers  around  the  interface;  these  parameters  are  determined  by  imposing 
electrostatic  consistency  among  the  perturbations,  Vh  and  the  charges,  n/( 
induced  at  each  layer  [10], 

A  crucial  point  coming  out  of  this  analysis  is  that  small  transfers  of  charge 
between  the  two  crystals  induce  large  electrostatic  potentials  [10].  Based  on 
this  fact,  we  can  use  a  first-order  approximation  and  assume  that  the  charges 
of  both  crystals  are  unaltered  except  for  two  layers  at  the  IF  (see  fig.  1),  in 
such  a  way  that  the  electrostatic  potentials  of  both  crystals  are  lined  up. 

Fig.  1  describes  this  first  approximation.  In  the  bulk  of  an  ionic  semicon¬ 
ductor  there  is  a  non-zero  charge,  Q,  on  each  anion  and  -Q  on  each  cation 
creating  an  electrostatic  potential.  The  alignment  of  the  mean  electrostatic 
potentials  of  semiconductors  a  and  fi  is  obtained  when  the  charges  on  the 
two  IF  layers  are  -(3 Qa  +  Qp)l 4  and  (Qa  +  3 Qf)/4  as  fig.  I  shows  (distance 
between  layers  are  assumed  to  be  the  same  in  all  cases).  This  is  for  two 
semiconductors  with  no  common  anion  or  cation  as  is  the  case  of  InAs-GaSb. 
For  GaAs-AlAs,  where  there  is  a  common  ion,  the  sequence  of  charges  is 

+  Qd)£ot)Qt, ...  • 

Tn  a  first  approximation  we  determine  V,. . .  V5  in  order  to  have  the  sequence 
of  charges  just  discussed  (the  IF  is  assumed  not  to  perturb  deeper  layers). 

The  computation  of  local  charge  densities  requires  the  calculation  of  the 
local  densities  of  states  on  each  layer  parallel  to  the  IF.  For  this  purpose  we 
compute  the  Green  function  of  the  system  projected  on  each  layer  by  means 
of  a  tight-binding  approach  where  the  spin-orbit  interaction  is  included  in  a 


Fig.  I.  Schematic  representation  of  charges  and  electrostatic  potentials  at  the  IF  in  the  case  of 
a  hctcroiunction  with  no  common  anion  or  cation. 
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sp’s*  basis  (details  of  this  calculation  will  be  published  elsewhere). 

The  relative  position  of  the  two  semiconductor  is  given  usually  interms  of 
the  offset  d£^0)  between  the  tops  of  the  valence  bands.  If  the  potential,  V$0' 
-  V^0),  we  have  calculated  in  the  first  approximation  coincides  with  the  differ¬ 
ence  between  the  mean  potentials  of  the  two  semiconductors,  having 
excluded  the  Hartree  potentials,  no  net  flux  of  charge  would  take  place  and 
the  solution  calculated  in  the  first  approximation  appears  to  be  already  con¬ 
sistent.  However,  this  is  not  usually  the  case  and  a  full  self-consistent  calcu¬ 
lation  is  necessary;  changes  in  V}  -  V,  would  give  the  change  in  AE$,0)  due 
to  the  self-consistent  calculation. 


3.  Results 

We  have  applied  the  method  described  in  section  2  to  the  heterojunction 
GaAs-AIAs,  InAs-GaSb  and  GaSb-InAs  along  the  (100)  direction.  Techni¬ 
cal  points  have  been  treated  as  in  ref.  (10]  except  for  the  integration  in  k/ 
that  we  have  performed  by  means  of  16  special  points  in  the  irreducible  part 
of  the  two-dimensional  Brillouin  zone.  Fig.  2  shows  the  local  density  of  states 
on  several  atomic  layers  both  around  the  IF  and  at  the  bulk  for  the  GaSb- 
inAs  heterojunctions  as  an  example  of  the  type  of  results  we  get.  Numerical 
results  for  the  potentials  at  the  IF  and  the  band  offset  are  given  in  table  1 
both  for  the  first  step  of  our  procedure  and  for  the  final  consistent  results.  It 
is  of  interest  to  note  that  the  final  consistent  calculation  is  rather  close  to  the 


Fig.  2.  Local  density  of  states  for  bulk  (8)  and  interface  (IF)  layers  as  a  function  of  the  energy 
Ain  eV)  in  the  GaSb-InAs  heterojunction. 


C.  Plaiero  el  al.  /  Electronic  structure  of  (100)  semiconductor  heterojunctkms 

Table  1 

Atomic  levels  around  the  IF  V,  and  valence  band  offset  AE,  (in  eV)  for  three  different 
heterojunctions;  the  magnitudes  with  superscript  have  been  obtained  in  the  first  step  of  the 
calculation  (see  text);  the  magnitudes  without  superscript  ate  the  final  consistent  values;  Vf" 
and  V,  are  taken  zero 


As  Ga  As  A!  As 

AstnSbGaSb 

SbGaAsInAs 

Vt" 

-0.011 

-0.086 

-0.343 

Vs 

0.025 

0.024 

-0.480 

V)"' 

-0.185 

-0.261 

-0.449 

v y 

-0.219 

-0.299 

-0.383 

vl"' 

-0.308 

-0.187 

-0.450 

vt 

-0.286 

-0.201 

-0.431 

vp 

-0.251 

-0.170 

-0.335 

Vy 

-0.248 

-0.044 

-0.341 

AEp 

-0.505 

0.610 

-1.115 

AE, 

-0  502 

0.736 

-1.121 

results  given  by  the  first  step:  similar  conclusions  have  been  reached  in  the 
metalfsemiconductor  junction  [10],  As  stated  above,  this  is  due  to  the  small 
transfer  of  charge  appearing  between  both  crystals. 

Comparing  now  our  results  with  the  experimental  information  we  find; 

(i)  GaAs-AlAs.  In  this  case  A £v  has  been  accepted  long  since  to  be  0.25 
eV,  a  result  obtained  by  extrapolating  [14]  for  x  -*  1  the  experimental  data 
for  GaAs-Ga,  _  ,AI,As.  However,  recent  experimental  information  [11) 
yields  d£v  *=  0.65  eV  and  a  heav)  hole  effective  mass  mj  *  0.34.  Since  the 
tight-binding  parameters  we  use  in  our  calculation  must  give  good  energy 
levels  at  high  symmetry  points,  we  cannot  get  any  desired  value  of  mj.  The 
best  value  we  get  is  mj  *  0.41;  then,  we  obtain  the  result  2l£v  =  0.50  eV 
given  in  table  1.  The  comparison  with  experimental  results  requires  an  estima¬ 
tion  of  the  error  involved  in  our  approximations.  We  think  that  this  uncer¬ 
tainty  must  be  of  the  order  of  0. 1  eV  which  would  imply  satisfactory  agree¬ 
ment  between  our  result  and  the  most  recent  experimental  value. 

(ii)  InAs-GaSb  In  this  case  the  most  interesting  result  is  that  we  find  a 
difference  of  0.4  eV  between  the  band  offset  of  InAs-GaSb  and  GaSb-InAs. 
The  difference  is  high  enough  to  allow  for  the  experimental  determination  of 
either  the  particular  IF  appearing  in  actual  samples,  or  the  electric  field 
associated  with  two  different  IFs  appearing  at  a  superlattice. 
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BAND-OFFSET  CALCULATIONS  OF  ( 10 0 ) -SEMICONDUCTOR 
HETEROJUNCTIONS:  A  SIMPLE  CONSISTENT  METHOD 

G.Platero,  J. SSnchez-Dehesa ,  C.Tejedor  y  F. Flores 
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Spain 


A  simple  consistent  tight-binding  method  is  presented 
to  analyse  the  electronic  structure  and  the  band-offset  of 
ionic  semiconductor  interfaces.  Our  model  is  an  extension  of 
the  model  based  on  the  charge  neutrality  level  concept,  and  our 
procedure  allows  us  to  calculate  the  induced  interface  dipole 
between  both  semiconductors.  We  present  results  for  GaAs-AlAs, 
GaSb-InAs  and  InAs-GaSb  (100)  hetero junctions .  Our  results 
yield  a  reasonable  agreement  with  the  experimental  evidence, 
and  show  different  band-offsets  for  the  GaSb-InAs  and  InAs-GaSb 


hetero junctions . 


1.  INTRODUCTION 


Many  electronic  properties  of  superlattices  depend 

crucially  on  the  band  offsets  between  the  different  semiconductors 

forming  the  junction.  Much  experimental  and  theoretical  work 

has  been  addressed  to  measuring  and  calculating  the  band  dis- 

1  2 

continuities  at  the  heterojunctions  '  . 

The  theoretical  work  can  be  classified  into  two  main 

groups.  In  one  class  of  papers,  people  have  tried  to  qive  simple 

prescriptions  to  calculate  the  band  discontinuities  ^  using 

bulk  properties  of  each  semiconductor.  One  of  the  most  promising 

ideas  in  this  kind  of  approach  was  introduced  some  time  ago  in 

Ref. 5,  and  recently  reformulated  and  discussed  by  Tersoff  ®  and 
2 

Margaritondo  .  The  main  concept  behind  this  model  is  the  charge 

neutrality  level,  the  energy  up  to  which  electronic  levels  are 

filled  in  a  free  surface.  Band  alignments  between  semiconductors 

are  obtained  by  aligning  the  two  charge  neutrality  levels  of 

both  semiconductors:  the  idea  supporting  this  "rule  of  thumb"  is 

that,  otherwise,  a  strong  electrostatic  dipole  would  be  induced 

at  the  interface  trying  to  restore  the  equilibrium  position 

defined  by  the  charge  neutrality  levels  alignment.  The  recent 

2 

results  presented  by  Margaritondo  give  a  strong  support  to 

this  kind  of  approach,  although  the  agreement  between  theory 

and  experiment  is  by  no  means  perfect. 

In  a  second  class  of  papers,  some  people  have  tried  to 

analyse  heterojunctions  by  means  of  a  full  self-consistent 
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calculation  using  a  pseudopotential  approach  .  No  doubt  this 
is  the  more  exact  approach,  and  were  the  calculations  simpler 
to  perform  and  less  computer  time-consuming,  we  would  have  in 


this  method  the  more  appropriate  theoretical  recipee  to  analyse 
hetero j  unctions . 

The  aim  of  this  paper  is  to  present  a  method  to  calculate 
hetero junctions  filling  the  gap  between  the  two  kind  of  work 
referred  to  above.  Here  we  present  a  theoretical  model  to 
calculate  the  hetero junction ,  based  on  a  tight-binding  method 
supplemented  with  an  appropriate  consistency;  this  method  offers 
a  step  forward  compared  with  the  simpler  models  discussed  above. 
Moreover,  and  this  is  the  crucial  idea  behind  our  approach,  our 
model  is  a.i  extension  of  the  model  based  on  the  charge  neutrality 
level  concept.  In  other  words,  the  consistency  of  our  calculation 
tries  to  take  into  account  the  charge  neutrality  conditions  and 
the  induced  interface  dipole  that  gives  the  basis  to  the  simple 
model  referred  to  above.  Let  us  mention  that  similar  ideas  have 
been  put  forward  for  the  metal-semiconductor  junctions  10  12 .  In 
this  case,  a  tight-binding  model  supplemented  with  a  Hartree 
consistency  has  been  recently  used  to  calculate  the  Schottky 
barrier  heights.  Results  have  been  encouraging,  showing  a  good 
agreement  with  the  experimental  evidence  and  other  more  elaborate 
approaches . 

In  §2  we  discuss  our  model,  our  method  of  calculation 
and  the  consistency.  In  §3  we  present  and  discuss  our  results 
as  applied  to  GaAs-AlAs,  GaSb-InAs  and  InAs-GaSb  (100)  hetero¬ 
junctions. 

2.  MODEL,  METHOD  OF  CALCULATION  AND  CONSISTENCY 

The  electronic  band  structures  of  the  different  semi¬ 
conductors  are  described  by  means  of  a  tight-binding  model,  using 


sp  s*  hybrids  in  each  atom  and  including  interaction  para¬ 
meters  that  extend  up  to  first  neighbours.  Table  1  gives  the 
parameters  used  in  our  calculation  for  GaAs,  AlAs,  GaSb  and  InAs. 
In  our  model,  we  assume  that  the  two  crystals  match  perfectly 
and  that  the  bulk  parameters  extend  up  to  the  interface,  except 
for  some  diagonal  perturbations  introduced  at  some  interface 
layers  to  impose  consistency  (see  below).  In  our  calculation  we 
have  also  included  the  spin-orbit  interaction  that  appears  to 

be  important  for  heavy  atoms  like  Sb.  The  parameters  given  in 

3 

Table  1  are  similar  to  the  ones  proposed  by  Vogl  et  al.  ,  except 
by  small  corrections  introduced  to  allow  for  the  spin-orbit 
interaction  (the  parameters  given  in  Table  1  have  been  chosen 
to  give  the  right  levels  at  specific  symmetric  points).  The 
parameters  defining  the  semiconductor-semiconductor  interactions 
are  associated  with  the  last  layer  of  the  crystals:  they  are 
defined  by  taking  the  interaction  parameters  for  the  crystal 
formed  by  the  corresponding  anion  and  cation. 

Having  defined  our  interactions  (except  for  the  diagonal 
perturbations  near  the  interface),  we  have  a  well-known  Hamilto¬ 
nian.  Our  next  task  is  to  discuss  how  we  calculate  the  interface 
properties:  charge,  density  of  states,  etc. 

In  our  method  of  calculation  we  follow  the  procedure 
explained  in  Ref. 14.  In  few  words,  our  method  consists  in  looking 
for  the  surface  components  of  the  Green  function  of  the  system. 

By  means  of  a  decimation  technique,  we  can  project  the  bulk 
components  of  the  Green  function  onto  a  few  layers  of  the  inter¬ 
face.  In  our  actual  case,  we  have  projected  all  the  electronic 
bulk  structure  of  each  semiconductor  onto  the  last  three  layers 


of  the  same  semiconductor,  the  reason  being  that  calculations 
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in  free  surfaces  show  that  only  the  last  three  layers  of  the 
semiconductor  are  perturbed  by  the  surface.  Thus,  for  GaSb-InAs, 
this  means  that  the  whole  system  is  projected  onto  six  layers, 
three  for  GaSb  and  three  for  InAs.  For  GaAs-AlAs,  the  same 
procedure  yields  only  five  layers,  since  there  is  a  central  As- 
atom  belonging  to  the  two  semiconductors  at  the  same  time.  This 
means  that  the  whole  Hamiltonian  is  reduced  to  an  effective 
matrix-Hamiltonian  having  60x60  elements  in  GaSb-InAs  and  50x50 
elements  in  GaAs-AlAs  (including  spin  and  five  orbitals  per 
layer).  From  this  matrix,  we  can  analyse  the  different  properties 
of  the  interface,  once  we  have  introduced  at  each  layer  the 
diagonal  perturbations  to  be  related  to  the  interface  consistency. 

Let  us  consider  now  how  to  get  consistency  in  our  problem. 
We  present  our  discussion  in  two  steps,  starting  with  the  InAs- 
GaSb  heterojunction  (GaAs-AlAs  will  be  discussed  later  on) . 

In  Figure  1  we  show  an  ideal  InAs-GaSb  hetero junction 
after  averaging  the  charges  in  the  sense  parallel  to  the  surface. 
For  InAs  we  assume  that  the  charges  per  atom  inside  the  crystal 
are  3-B(In)  and  5+B(As),  while  for  GaSb  the  charges  are  assumed 
to  be  3-«(Ga)  and  5+a(Sb).  The  ideal  case  shown  in  Figure  1  is 
defined  by  assuming  that  the  electrostatic  potential  inside  each 
crystal  is  the  same  as  in  the  bulk,  and  that  both  mean  electros¬ 
tatic  potentials  are  equal.  Notice  that  for  these  conditions, 
the  charges  in  each  layer  coincide  with  the  bulk  charges  except 
for  the  last  layers  of  each  semiconductor  (in  Figure  1,  As  for 
InAs  and  Ga  for  GaSb) .  Charges  in  those  layers  must  be  modified 
in  order  to  have  both  mean  electrostatic  crystal  potentials  at 


the  same  level.  This  condition  yields  (see  Figure  1) : 
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The  ideal  situation  shown  in  Figure  1  allows  us  to  define  a 
starting  point  for  which  the  different  parameters  appearing  in 
the  Hamiltonian  are  given  by  the  bulk  parameters  of  each 
crystal  the  two  mean  levels  of  both  crystals  are  corrected 
in  order  to  take  into  account  the  condition  of  equality  for  the 
mean  electrostatic  potentials.  That  means  that  for  the  charges 
shown  in  Figure  1,  we  can  expect  to  have  at  each  layer  the  ideal 
parameters  defined  for  a  bulk  Hamiltonian:  thus,  assume  we  take 
the  parameters  of  Table  1  for  InAs;  then,  for  GaSb  we  must  take 
for  the  atomic  levels  of  each  atom  the  values  of  Table  1  but 
corrected  by  a  constant.  A,  putting  the  mean  electrostatic 
potential  of  both  crystals  at  the  same  level.  This  value,  A, 
can  be  taken  from  an  independent  information  about  the  surfaces 
of  each  crystal  on  the  other  hand,  the  hopping  inter¬ 

action  between  the  orbitals  of  Ga  and  As,  at  the  interface,  are 
taken  from  the  values  corresponding  to  GaAs. 

Once  we  have  defined  this  initial  situation,  we  now 
allow  for  changes  in  the  charges  and  potentials  of  the  layers 
around  the  interface.  Notice  that  for  the  initial  potentials 
defined  above,  we  can  obtain  the  charges  at  each  layer  by 
solving  the  total  Hamiltonian:  in  general,  the  charges  at  each 
layer  do  not  coincide  with  the  initial  charges  shown  in  Figure 


1.  In  our  consistent  procedure,  we  allow  for  changes  in  the 
potentials  of  layers  2,  3,  4,  5  and  6  (taking  layer  1  as  a 
reference  layer),  and  assume  that  the  mean  level  of  GaSb  is 
shifted  by  the  same  value  of  layer  6.  Then,  the  changes  in  the 
potentials,  6V2,  6V3,  ...,  6Vg  are  related  to  the  charges  of 
each  layer  through  a  kind  of  consistent  equation.  This  means  that 
are  defined  as  the  electrostatic  potential  created  by  the 
charges  6n^,  as  measured  with  respect  to  the  ideal  case  shown 
in  Figure  1.  Thus, 

<5V2  =  ad6n^ 

<SV3  =  ad6n^  +  ad(<5n^+6n2) 

(2) 

SVg  =  ad6njL  +  ad(6n1+6n2)  +  ...  +  ad  (  6n3  +  6n2+6n3  +  6n4+6n5)  , 

where  d  is  the  interlayer  distance  (assumed  to  be  a  constant) 
and  a  =  4tt/a,  A  being  the  area  per  surface  atom. 

Notice  that  in  our  actual  case  we  have  found  that 

equation 

<5n^  +  6n2  +  ...  +  6n^  +  6ng  =  0  (3) 

is  verified  with  a  very  good  approximation  (this  means  that  the 
interface  perturbation  only  penetrates  three  layers  into  each 
semiconductor) . 

Eqs.(2)  define  our  consistent  conditions,  and  allow  us 


to  obtain  the  final  perturbations,  6V^,  and  the  final  charges 
at  the  interface. 

The  GaAs-AlAs  interface  can  be  analysed  in  similar 
terms.  Figure  2  defines  the  starting  point  for  this  hetero¬ 
junction.  Here,  we  assume  that  the  electrostatic  crystal  potential 
for  each  crystal  is  not  modified  up  to  the  central  layer  of  As. 

In  this  initial  situation,  we  assume  both  crystals  to  have  the 
same  mean  electrostatic  potential  (as  in  the  previous  case, 
that  means  that  the  atomic  levels  of  one  crystal,  say  GaAs,  have 
to  be  shifted  by  a  constant.  A,  with  respect  to  the  atomic 
levels  of  AlAsJ.  The  greater  difference  with  the  previous  case 
appears  for  the  central  atom:  as  Figure  2  shows,  the  ideal 
potential  associated  with  each  ideal  crystal  yields  a  different 
potential  for  the  central  As-atom.  We  define  the  starting  case 
by  taking  for  this  central  atom  the  averaging  of  the  two 
potentials  associated  with  each  crystal.  Then,  the  starting 
potential  looks  like  the  one  shown  by  the  continuous  line  of 
Figure  2;  for  this  potential  we  have  to  introduce  some  transfer 
of  charge  among  layers  2,  3  and  4.  It  is  easy  to  see  that  the 
initial  charges  defined  by  the  initial  electrostatic  potential 
are  given  by  (see  Figure  2): 


Once  we  have  defined  the  initial  case,  we  look  for 


consistency  allowing  for  changes  in  the  potentials  at  layers 
2,  3,  4  and  5.  Then,  V2...  and  V,.  are  determined  consistently 
by  writing  the  following  Hartree  equations: 

6V2  =  adfin.^ 


<5V3  =  ad  6n^  +  ad(6n1+6n2) 


(5) 


6V5  =  adSn^j^  +  ad(6n1+6n2)  +  ...  +  ad  ( Sn^t. .  ,  +  fin^ )  , 

where  d  is  the  interlayer  distance,  and  a  =  4n/A,  A  being  again 
the  area  per  surface  atom. 

Eqs.  (5)  define  the  consistent  equations  for  the  GaAs- 
AlAs  interface.  Notice  that  6n^  must  be  consistent  with  the 
values  calculated  by  introducing  6Vi  in  the  total  Hamiltonian. 

3.  RESULTS  AND  DISCUSSION 

We  have  analysed  the  ideal  interfaces  associated  with 
the  following  heterojunctions :  GaAs-AlAs,  InAs-GaSb  and  GaSb-InAs. 
(for  InAs-GaSb  the  interface  corresponds  to  atoms  of  As  and  Ga, 
while  for  GaSb-InAs,  it  corresponds  to  Sb  and  In).  Technical 
points  have  been  treated  as  in  Ref. 10,  except  for  the  integration 
in  the  momentum  parallel  to  the  interface,  which  we  have 
performed  by  means  of  16  special  points  in  the  irreducible  part 
of  the  two-dimensional  Brillouin  zone. 
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In  our  calculations,  we  have  used  the  following  values 
16,17. 

A  (GaAs-AlAs )  =  0.26  eV 
A  (GaSb-InAs)  =0.78  eV. 

The  different  diagonal  perturbations  calculated  consistently 
for  the  different  interfaces  are  given  in  Table  2.  For  GaAs-AlAs 
we  have  obtained  an  interface  dipole  of  -0.13  eV  (the  value  of 
V5  in  Table  2),  and  the  following  discontinuity  between  the 
top  of  the  valence  bands: 

AEy  =  A  +  v5  =  0.41  eV  (GaAs-AlAs)  .  (6) 

(A  is  the  valence  band-offset  for  no  interface  dipole) . 

Figure  3  shows  the  local  density  of  states  on  several 
atomic  layers  around  the  interface  of  a  consistent  GaAs-AlAs 
junction.  For  GaSb-InAs,  we  have  calculated  an  interface  dipole 
of  -0.26  eV  (the  value  of  Vg  in  Table  2),  and  the  following 
discontinuity  in  the  top  of  the  valence  bands  21 : 

AEV  =  A  +  Vg  =  0.52  eV  (GaSb-InAs)  ,  (7) 

while  for  InAs-GaSb,  our  calculations  yield  for  the  interface 
dipole  0.53  eV,  and 

E v  =  A  +ve  *  1 • 3 1  eV  (InAs-GaSb)  .  (8) 


Figures  4  and  5  show  the  local  density  of  states  on 
several  layers  around  the  interfaces  of  GaSb-InAs  and  InAs-GaSb, 
respectively. 

The  important  difference  between  the  GaSb-InAs  and 
InAs-GaSb  interfaces  comes  from  the  different  interface  dipoles 
appearing  in  both  cases  (-0.26  eV  and  0.53  eV,  respectively). 

Note  that  we  can  go  from  one  to  the  other  interface  by 
substituting  an  In  for  a  Ga  atom.  The  different  dipoles  appearing 
for  both  interfaces  are  associated  with  the  sequences 
. . .-Sb-Ga-As-. . .  and  . . .-Sb-In-As-. . .  of  both  junctions.  Thus, 
the  different  atoms  of  Ga  and  In  are  controlling  the  way  in  which 
the  electronic  charge  is  transferred  from  Sb  to  As  and,  finally, 
from  one  to  other  crystal.  The  crucial  point  is  that  the  InAs 
bond  is  associated  with  an  important  transfer  of  electronic 
charge  towards  As;  when  In  is  substituted  for  Ga,  part  of  this 
charge  is  transferred  from  As  to  sb  creating  an  interface  dipole 
decreasing  the  valence  band  discontinuities. 

Comparing  now  our  results  with  other  experimental  and 
theoretical  information  we  find; 

(i)  GaAs-AlAs.  For  this  case,  AEy  has  been  accepted  long 

18 

since  to  be  0.25  eV  ,  a  result  obtained  by  extrapolating  for 
x  -*■  1,  the  experimental  data  for  GaAs-Ga,  A1  As.  However, 

i.  A 

19 

recent  experimental  information  yields  AEy  =  0.65  eV.  Our 
result  of  0.41  eV  is  in  a  reasonable  agreement  with  that 
information;  note  in  that  respect  that  the  accuracy  of  our 
calculation  is  around  ±0.05  eV,  and  that  the  error  involved  in 
the  approximations  of  our  method  is  estimated  to  be  around  ±0.1  eV 

(ii)  InAs-GaSb  and  GaSb-InAs.  In  this  case  we  find  the 


following  band-offset  discontinuities:  1.31  eV  (InAs-GaSb)  and 

0.52  eV  (GaSb-InAs).  These  values  have  to  be  compared  with  the 

20 

experimental  data  of  0.57  eV  .  According  with  the  errors 
involved  in  our  method  (see  above) ,  we  conclude  that  only  the 
value  found  for  the  GaSb-InAs  interface  seems  to  be  in  good 
agreement  with  the  experimental  evidence.  This  seems  to  suggest 
that  the  actual  interface  of  GaSb  and  InAs  correspond  to  the 
GaSb-InAs  case.  However,  a  word  of  caution  must  be  put  here, 
since  the  calculations  presented  in  this  paper  have  been 
performed  neglecting  the  mismatch  between  both  crystals;  the 
effect  of  this  mismatch  ought  to  be  analysed  before  definitive 
conclusions  are  drawn. 

Let  us  finally  comment  that  the  different  band-offset 
discontinuities  found  for  the  GaSb-InAs  and  InAs-GaSb  interfaces 

o 

seem  to  be  at  variance  with  results  reported  by  Picket  et  al.  . 
In  this  sense,  it  is  of  interest  to  comment  that  in  this  last 
calculation  the  number  of  layers  used  in  each  superlayer  of  each 
semiconductor  is  small  enough  (only  5)  to  allow  for  an  important 
interaction  between  different  interfaces;  note  in  this  respect 
that  the  interface  perturbation  penetrates  at  least  three  layers 
in  each  semiconductor. 
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Figure  Captions 

Figure  1.  Electrostatic  potential  of  an  ideal  InAs-GaSb  hetero¬ 
junction  after  averaging  the  charges  in  the  sense 
parallel  to  the  interface. 

Figure  2.  As  Fig.l  for  a  GaAs-AlAs  hetero junction. 

Figure  3.  Local  density  of  states  on  different  atomic  layers 
around  the  interface  of  a  GaAs-AlAs  junction. 


Figure  4.  As  Fig. 3  for  GaSb-InAs. 
Figure  5.  As  Fig. 3  for  InAs-GaSb. 
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The  initial  stages  of  the  Schottky-barrier  formation  for  an  Al— Si  contact  are  theoretically 
analysed.  Four  different  geometrical  configurations  are  discussed,  and  their  charge  neutrality 
level,  interface  Fermi  level  and  local  density  of  states  are  calculated.  We  conclude  that  the 
Schottky  barrier  is  formed  with  a  monolayer  and  that  the  charge  neutrality  level  at  the  junction 
is  a  function  of  the  metal-semiconductor  interaction. 


The  mechanism  of  the  Schottky-barrier  formation  is  a  matter  of  current 
controversy  [1],  Several  models  have  been  proposed  in  the  last  few  years. 
The  defect  model  [2]  seems  to  be  very  appropriate  for  ionic  semiconductors; 
the  problem  with  it  is  related  to  the  number  of  defects  necessary  to  pin  the 
Fermi  level  at  the  interface.  A  different  mechanism  based  on  the  concept  of 
the  charge  neutrality  level  at  the  interface  was  first  proposed  by  Tejedor  et 
al.  [3];  recently  Guinea  et  al.  [4]  and  Tersoff  (5]  have  extended  this  approach, 
which  seems  to  be  particularly  appropriate  for  abrupt  interfaces  with  a  well- 
defined  separation  between  the  metal  and  the  semiconductor.  In  the  mecha¬ 
nism  of  this  model,  charge  neutrality  conditions  fix  the  interface  Fermi  level 
at  an  energy  quite  close  to  the  Fermi  level  for  a  free  surface.  Detailed  calcu¬ 
lations  using  a  self-consistent  tight-binding  method  [4,  6]  have  shown  that  the 
Fermi  level  of  the  free  surface  cah  be  slightly  shifted  in  a  junction,  depending 
on  the  interaction  between  the  last  layers  of  the  metal  and  the  semiconductor. 
Experimental  results  tend  to  support  these  findings;  thus,  Margaritondo  et 
al.  [7]  have  found  that  for  Al,  Ga  and  In  on  jSi,  the  Fermi  level  at  the  interface 
moves  downwards,  with  respect  to  the  free  surface,  by  0.2,  0.3  and  0.4  eV, 
respectively. 

In  this  paper  we  consider  the  initial  stages  in  the  Schottky  barrier  forma¬ 
tion;  in  particular,  we  have  analysed  the  case  of  a  monolayer  of  Al  on  Si. 
This  case  has  been  studied  experimentally  [8],  •  id  the  evidence  seems  to 
indicate  that  the  Schottky  barrier  is  already  formed  with  only  a  monolayer. 
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In  previous  papers,  Chang  and  Schluter  [9],  and  Chelikowsky  [10]  have 
analysed  theoretically  the  same  problem  for  a  monolayer  in  a  Lander  substi¬ 
tutional  geometry  and  in  other  surface  configurations.  One  of  these  latter 
configurations  corresponds  to  an  Al  atom  on  top  of  a  Si  surface  atom,  another 
corresponds  to  an  Al  atom  placed  on  a  three-fold  coordinated  site.  The  ana¬ 
lysis  of  Chang  and  Schluter  [9]  tends  to  support  a  three-fold  coordinated  site 
for  the  Al  atom.  However,  the  theoretical  results  of  these  authors  show  that 
the  Fermi  level  at  the  surface  is  shifted  to  higher  energies  with  respect  to  the 
Fermi  level  of  the  free  surface  by  0.2  eV,  this  result  being  in  conflict  with 
the  fact  that,  for  a  full  monolayer,  the  batrier  is  completely  formed  and  that 
the  Fermi  level  for  the  junction  is  located  at  lower  energies  than  the  Fermi 
level  for  the  free  semiconductor  surface. 

In  this  work,  we  reconsider  the  problem  of  a  -monolayer  of  AJ  on  Si,  and 
discuss  four  different  configurations:  the  top  position,  two  three-fold  coordi¬ 
nated  positions  which,  following  Northrup  [11],  will  be  called  the  T3  and  H4 
configurations  (Northrup  has  shown  that  the  H4  site,  with  an  Al  atom  placed 
in  a  three-fold  site  having  a  Si  atom  underneath  in  the  second  layer,  is  the 
most  stable  for  a  VI  x  y/Z  configuration),  and  a  Lander  substitutional  model. 
In  our  calculations,  we  have  followed  ref.  [4],  and  have  calculated  the  electro¬ 
nic  density  of  states  by  means  of  a  self-consistent  tight-binding  procedure.  In 
table  1  we  give  the  parameters  used  to  calculate  the  Si  band  structure  (these 
parameters  are  similar  to  the  ones  introduced  by  Pandey  and  Phillips  [12]), 
the  AI-AI  interactions  (a  first-neighbours  interaction  calculated  with  an  in¬ 
verse  square  law)  and  the  Si-Al  interactions  (calculated  by  taking  a  geometri¬ 
cal  average  between,  the  interactions  for  the  nearest  neighbours  of  Si-Si  and 
AI-AI). 

Figs.  1  and  2  show  the  local  density  of  states  projected  on  the  Al  layer  as 
calculated  with  our  procedure  for  the  top,  Hj,  T4  and  the  substitutional  con¬ 
figurations.  For  completeness,  we  have  calculated  results  for  a  clean  surface 
(in  this  latter  case,  the  density  of  states  is  projected  on  the  last  Si  layer).  Our 
results  show  that:  (i)  in  the  top  position,  the  local  density  of  states  in  Al  is 
more  localised  near  the  Fermi  level,  corresponding  to  a  weaker  interaction 
between  AJ  and  Si;  (ii)  the  substitutional  case  is  similar  to  the  clean  Si  surface, 
but  presenting  empty  surface  states;  (iii)  the  Hj  and  T4  positions  have  similar 
densities  of  states  corresponding  to  similar  interactions  between  Al  and  Si. 

These  results  are  in  correspondence  with  the  position  of  the  Fermi  level 
at  the  interface.  Table  2  gives  the  shift  of  the  Fermi  level  for  the  different 
interfaces  as  measured  from  its  position  at  the  free  surface.  In  table  2  we  also 
show  the  Fermi  level  for  the  case  in  which  no  transfer  of  charge  is  allowed 
between  the  metal  and  the  semiconductor  (this  latter  case  defines  the  charge 
neutrality  level  for  each  interface).  Our  results  show  that  the  charge  neutrality 
level  of  each  interface  is  different  and  dependent  on  the  interaction  between 
the  metal  and  the  semiconductor.  In  the  top  position,  the  charge  neutrality 
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Table  1 

Different  parameters  defining  the  Hamiltonian  used  in  the  calculation  (the  Si-AI  distance  is 
taken  the  same  in  all  the  cases) 


Si-Si 

Si-AI 

Al-Al 

E, 

-3.30 

-0.65 

Ep 

1.10 

- 

2.65 

(sso), 

-2.08 

-1.30 

-0.32 

(spo)i 

1.77 

2.42 

1.20 

(ppo)l 

2.52 

2.99 

1.28 

(PPtr)i 

-0.33 

-0.16 

0.0 

(PP°h 

0.58 

- 

“ 

(PP*)j 

-0.10 

- 

ENERGY  (eV) 

Fig.  I.  Local  density  of  states  on  the  last  Si  layer  for  (a)  the  top  position,  (b)  a  Lander  substitu¬ 
tional.  and  (c)  a  clean  S  surface.  Er  a  Fermi  level.  The  arrow  in  (a)  shows  the  specific  surface 
peak  found  in  photoemission  experiments  (7). 
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level  is  shifted  downwards  in  energy,  while  for  larger  interactions  between 
the  metal  and  the  semiconductor  (H3  and  T4  positions)  that  level  is  shifted 
upwards.  Let  us  remark  that  the  Fermi  level  does  not  exactly  coincide  with 
the  charge  neutrality  level:  the  main  reason  is  that  for  the  full  solution  the 
charge  neutrality  case  has  to  be  relaxed,  allowing  for  some  transfer  of  charge 
between  the  metal  layer  and  the  last  Si  layer. 

Let  us  comment  at  this  point  that  the  absolute  values  we  have  calculated 
for  the  Fermi  level  shift  at  different  interfaces  are  too  large  due  to  the  Pan- 


IV* 


-12  -8-4  0 

ENERGY  lev  I 

Fig.  2.  A*  fig.  I  for  (a)  a  clean  surface,  (b)  a  H3  position,  and  (c)  a  T,  position. 


Table  2 

The  Fermi  level,  the  charge  neutrality  level  and  the  Fermi  level  shift  for  different  geometries 
(the  aero-energy  level  it  located  at  the  top  of  the  Si  valence  band) 


Geometrical  position  Et  (eV) 


Si  (dean  surface)  0.S5 

Top  0.00 

Substitutional  0.25 

Hj  1.10 

T«  1.40 


Fermi  level  shift  (eV)  Charge  neutrality  level  (eV) 


W' 
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dey— Phillips  model  we  use  for  the  semiconductor.  As  discussed  in  ref.  [4], 
the  Pandey-Phillips  model  yields  values  of  the  Fermi  level  shift  which  are  too 
large  by  a  factor  of  two,  due  to  the  large  value  of  the  thermal  gap. 

It  is  also  interesting  to  note  that  the  results  calculated  here  for  the  Al  top 
position  are  in  good  agreement  with  the  ones  given  in  ref.  [4],  for  a  semi-in¬ 
finite  metal  having  the  atoms  of  the  last  layer  located  at  the  top  position,  too. 
This  shows  that  the  barrier  height  is  practically  formed  with  a  monolayer,  in 
agreement  with  experiments. 

Comparing  now  with  the  experimental  evidence  [8],  let  us  comment  that 
(i)  the  Fermi  level  is  shifted  to  lower  energies  for  a  well-formed  contact,  and 
that  (ii)  the  density  of  states  is  very  low  near  the  Fermi  level,  showing  a 
specific  peak  near  0.7  eV  below  the  Fermi  level.  These  results  are  in  reason¬ 
able  agreement  with  our  calculations  for  the  top  position.  In  this  case  the 
Fermi  level  is  shifted  by  0.55  eV  (a  value  too  large  by  a  factor  of  two),  and 
the  local  density  of  states  shows  a  similar  peak  to  the  one  found  in  the 
experimental  results  (see  fig.  la). 

We  conclude  that  the  Al  monolayer  is  completed  with  the  metal  atoms 
located  on  the  top  position.  Northrup's  results  suggest  that  a  V5  x  V5  Al 
monolayer  has  the  metal  atoms  in  a  T4  position.  These  two  conclusions  are 
not,  howeveT,  contradictory  if  we  consider  the  strong  interaction  operating 
between  Al  atoms  in  the  full  monolayer  case. 

Finally,  it  is  worth  stressing  our  results  for  the  charge  neutrality  level  of 
the  junction.  As  our  results  have  shown,  that  level  is  a  function  of  the  specific 
interaction  between  the  last  layers  of  the  metal  and  the  semiconductor.  This 
is  an  important  point  which  must  be  taken  into  account  to  explain  the  barrier 
heights  of  different  junctions  for  the  same  semiconductor. 
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A  simple  consistent  tight-binding  method  is  used 
to  analyse  the  initial  stages  of  the  Schottky  barrier 
formation  for  an  abrupt  Al-GaAs (100)  interface.  Two  different 
cases  are  discussed:  (i)  firstly,  we  consider  the  case 
of  an  A1  layer  with  the  Al  atoms  located  in  the  cation 

position  of  the  ideal  continued  GaAs(lOO)  structure;  (ii) 
in  the  second  case,  we  consider  a  full  Al  monolayer  with 

twice  the  number  of  Al  atoms,  having  a  geometry  similar 
to  the  (100) -face  of  an  Al  f.c.c.  structure.  Our  results 

show  that  the  Schottky  barrier  is  completely  formed  with 

the  full  monolayer,  yielding  a  barrier  height  of  0.55*0.05 
eV.  This  barrier  height  is  similar  to  the  one  calculated 
for  an  abrupt  roetal-GaAs ( 110 )  interface. 


1.  INTRODUCTION 


Different  models  of  Schottky  barrier  formation  have 

been  proposed,  and  it  is  clear  that  no  model  has  a  universal 

acceptance  since  no  single  model  can  explain  all  the  experimental 

data.  This  is  related  to  the  different  kind  of  interfaces 

appearing  in  a  semiconductor-metal  junction;  according  with  the 

1  2 

experimental  evidence  we  find:  (i)  abrupt  interfaces  '  with  a 

well  defined  separation  between  the  structures  of  the  metal 

and  the  semiconductor,  and  (ii)  reactive  interfaces ^  with 

the  metal  and  the  semiconductor  interdif fusing  and/or  forming  a 

new  chemical  compound  at  the  interface. 

The  different  models,  related  to  the  two  kinds  of 

interfaces  commented  on  above,  appear  to  continue  to  attract 

fairly  wide  support  in  the  scientific  community:  (i)  The  defect 
7-11 

model  ,  which  assumes  that  the  Schottky  barrier  formation 

results  from  the  creation  of  significant  number  of  defects  at 

or  near  the  interface,  and  (ii)  the  intrinsic  interface  states 
22-16 

model  ,  which  attributes  the  Schottky  barrier  height  to 

Fermi  level  pinning  by  the  metal  induced  density  of  states  near 
the  semiconductor  midgap.  Both  models  involve  states  at  the 
metal-semiconductor  interface  with  properties  that  are  almost 
independent  of  the  metal,  explaining  Schottky  barrier  heights 
that  are  virtually  independent  of  the  metals. 

There  is  a  general  agreement  in  the  fact  that  the 
mechanism  determining  the  Schottky  barrier  heights  depends  on 
the  kind  and  quality  of  the  interface.  Thus,  the  defect  model 
can  be  accepted  as  the  mechanism  controlling  the  Schottky 
barrier  height  li  the  number  of  defects  at  the  interface  is 


high  enough  (something  between  0.1%  and  10%  of  the  surface 

atoms,  a  subject  of  current  controversy). 

In  this  paper  we  analyse  the  initial  stages  of  the 

Schottky  barrier  formation  for  an  abrupt  Al-GaAs(lOO)  interface. 

The  quality  of  this  interface,  when  A1  is  epitaxially  grown  in 

GaAs(lOO),  appears  to  be  very  dependent  on  the  temperature  of 
17  18 

deposition  '  .  For  low  enough  temperature,  the  interface  is 

abrupt,  and  the  epitaxy  of  A1  is  observed  to  grow  with  different 

18 

orientations  on  GaAs  .  For  simplicity,  in  this  paper  we  only 

consider  the  case  for  which  Al(100)  growth  is  observed,  a  case 

that  would  be  expected  on  the  basis  of  simple  lattice  metals 

arguments  (see  below) .  According  with  the  previous  discussion, 

the  Schottky  barrier  formation  for  this  low  temperature  case 

can  be  expected  to  be  controlled  by  the  metal  induced  density  of 

states  near  the  semiconductor  midgap.  This  is  the  point  of  view 

taken  in  this  paper,  and  the  first  stages  of  the  Schottky 

barrier  formation  is  analysed  following  this  model. 

As  regards  the  different  theoretical  approaches  used  to 

analyse  the  ideal  abrupt  interfaces  and  the  metal  induced 

density  of  states  at  the  semiconductor  midgap,  let  us  mention: 

(i)  the  self-consistent  calculations  using  pseudopotential 

19 

theory  and  local  density  formalism  ,  and  (ii)  the  simple 

models  introducing  the  charge  neutrality  level  at  the  interface 
13  16 

'  .  In  this  last  context,  a  mechanism  of  Schottky  barrier 

formation  based  on  the  concept  of  charge  neutrality  level  was 
first  proposed  in  Ref. 13.  Here,  it  was  shown  that  for  a  covalent 
surface,  when  the  semiconductor  and  the  metal  are  coupled,  the 
semiconductor  density  of  surface  states  broadens  but  keeps  its 


centre  of  gravity  almost  constant.  Hence  the  interface's  Fermi 

level  almost  coincides  with  the  Fermi  level  for  the  free  surface 

For  ionic  surfaces,  like  the  (100)  of  GaAs,  things  are  more 

complicated:  in  this  case,  the  centre  of  gravity  of  the  surface 

bands  does  not  coincide  with  the  mean  level  of  the  optical 
20 

gap  ;  on  the  contrary,  due  to  the  net  charge  of  the  last 

atomic  layer,  the  centre  of  gravity  of  the  surface  bands  is 

21 

displaced  for  an  anion  surface  towards  the  valence  band 

However,  for  an  intimate  metal-semiconductor  junction,  the 

surface  band  of  the  ionic  surface  also  broadens  and  the  final 

density  of  interface  states  resembles  very  much  that  of  a 
21 

covalent  case  .  In  line  with  these  comments,  Tersoff  has 

22 

suggested  recently  that  the  Schottky  barrier  for  the  ionic 

face  of  a  semiconductor  can  be  also  obtained  by  using  the  same 

charge  neutrality  level  calculated  for  a  neutral  face  of  the 

same  semiconductor  (say,  a  (110) -surf ace) . 

In  this  paper  we  have  followed  an  alternative  approach, 

recently  used  succesfully  to  analyse  Schottky  barriers  for 

14  15 

covalent  semiconductors  '  .  This  is  a  tight-binding  method 

supplemented  with  an  additional  consistency  associated  with 

charge  neutrality  conditions.  The  idea  behind  this  approach  is 

13 

to  extend  the  charge  neutrality  level  concept  to  a  more 
sophisticated  calculation  which  takes  into  account  the  small 
variations  that  interface  dipoles  introduce  in  the  "rule  of 
thumb"  associated  with  the  charge  neutrality  approach. 

On  the  other  hand,  our  calculation  is  directed  to 
understanding  the  first  stages  of  the  Schottky  barrier  formation 


In  a  previous  work  ,  we  found  that  for  A1  on  Si  (111),  the 
barrier  height  is  completely  formed  for  a  monolayer  of  Al; 
this  suggests  that  a  similar  case  can  be  found  in  the  actual 
junction  we  discuss  here.  Our  final  results  show  that,  indeed 
for  a  complete  monolayer  of  Al  on  GaAs(lOO)  a  finite  density 
of  states  fill  the  semiconductor  gap,  suggesting  that  the 
barrier  is  completely  formed . 

In  §2  we  discuss  our  model,  our  method  of  calculation 
and  consistency,  and  in  §3  wer present  and  discuss  our  results 


2.  MODEL,  METHOD  OF  CALCULATION  AND  CONSISTENCY 


In  this  work  we  consider  the  problem  of  a  monolayer 
(and  a  half  of  a  monolayer)  of  Al  on  GaAs(lOO).  Figure  la 
shows  one  of  the  geometrical  configurations  we  have  analysed. 
Here,  one  of  the  Al  atoms  is  located  in  the  Ga  position  for 
the  ideal  continued  GaAs  structure.  For  this  position,  the  Al 
layer  forms  an  ideal  (100) -face  of  the  Al  metal  structure, 
except  for  other  Al  atoms,  located  in  the  centre  of  the  elemental 
face  formed  by  four  Al  atoms.  We  consider  two  structures:  (i) 
the  first  one  only  includes  Al  atoms  in  the  ideal  positions  of 
Ga,  (ii)  the  second  one,  the  full  monolayer  case,  has  twice 
the  number  of  Al  atoms,  and  is  similar  to  the  (100) -face  of 
the  Al  f.c.c.  structure.  For  this  case,  the  Al-(llO)  direction 
in  the  (lOO)-plane  is  parallel  to  the  GaAs- (010)  direction. 
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The  electronic  band  structure  of  GaAs  is  described  by 

3  .  24 

means  of  a  tight-binding  model,  using  sp  s  hybrids  in  each 

atom,  and  including  interaction  parameters  that  extend  up  to 

first  neighbours.  The  effect  of  the  interface  penetrates  up  to 

25 

the  third  layer  of  GaAs  ;  accordingly,  we  introduce  diagonal 
perturbations,  and  V2 ,  in  the  last  two  layers  of  the  semi¬ 
conductor.  These  perturbations  will  be  calculated  consistently 
as  discussed  below.  As  regards  Al,  we  introduce  interactions 

between  As  and  the  nearest  Al  atom  according  with  the  parameters 
24 

of  AlAs  .  The  Al-Al  interactions  are  the  same  used  in  Ref. 23, 
and  they  have  been  calculated  using  Harrison's  prescription  and 
an  inverse  square  law  for  the  dependence  on  the  Al-Al  distance. 
In  Table  I  we  give  the  different  parameters  defining  the  inter¬ 
face  Hamiltonian. 

In  our  method  of  calculation  we  follow  the  procedure 
developed  in  Ref. 14.  Summarising,  our  method  uses  a  Green 
function  approach  and,  by  means  of  a  decimation  technique, 
calculates  the  surface  components  of  the  Green  function  of  the 
system.  Thus,  we  can  project  the  bulk  components  of  the  whole 
Green  function  onto  a  few  layers  of  the  interface.  In  our 
actual  case,  we  have  projected  all  the  electronic  bulk  structure 
of  the  semiconductor  onto  the  last  three  layers,  say,  the  last 
As  layers,  and  the  next  Ga  and  As  layers.  In  this  way,  we 
reduce  the  whole  Hamiltonian  to  an  effective  one,  associated 
with  the  Al  layer,  and  the  three  last  layers  of  GaAs  (this 
yields  a  20x20  effective  matris  for  a  half  monolayer  of  Al , 
counting  five  orbitals  per  layer  and  a  25x25  matrix  for  the 


full  monolayer.  From  this  matrix,  we  can  analyse  the  different 
electronic  properties  of  the  interface,  once  we  introduce  at 
each  layer  (the  A1  layer  included)  the  diagonal  perturbations 
to  be  related  to  the  interface  consistency. 

Let  us  now  consider  how  to  get  the  interface  consistency 
in  our  problem.  In  Figure  2,  we  show  an  ideal  Al-GaAs  interface, 
after  averaging  the  charges  of  each  atom  in  the  sense  parallel 
to  the  surface.  For  GaAs  we  assume  that  the  charges  per  atom 
inside  the  crystal  are  3-a  (Ga)  and  5+a  (As) .  The  ideal  case 
shown  in  Figure  2  is  defined  by  assuming  that  the  electrostatic 
potential  near  the  interface  is  the  same  as  in  GaAs,  and  that 
the  electrostatic  potential  between  the  last  layers  of  A1  and 
As  coincides  with  the  one  for  an  AlAs  crystal;  this  means  that 
the  charge  per  atom  for  the  last  Al  layer  is  3-B/2,  while  the 
charge  per  atom  for  the  last  As  layer  is  5+a/2+B/2,  B  being  the 
transfer  of  charge  between  atoms  of  Al  and  As  in  an  AlAs  crystal 
Notice  that  this  ideal  case  allows  us  to  define  the  parameters 
for  Al  (as  given  in  Table  1)  by  referring  the  Al  levels  to  the 
As  ones  by  means  of  the  parameters  of  the  AlAs  crystal. 

The  ideal  situation  shown  in  Figure  1  allows  us  to 
define  a  starting  point  for  which  the  different  interactions 
appearing  in  the  Hamiltonian  are  given  by  the  parameters  of 
Table  1.  That  means  that  for  the  charges  shown  in  Figure  2,  we 
expect  to  have  the  parameters  given  in  Table  1. 

Once  we  define  this  initial  situation,  we  now  allow  for 
changes  in  the  charges  and  potentials,  VQ,  and  V2  (see 
Figure  2),  of  the  layers  around  the  interface.  Notice  that  for 


the  initial  potentials  defined  above  (for  this  case  we  define 


Vq  =  V1  =  V2  =  0)  we  can  obtain  the  charges  at  each  layer  by 
solving  the  total  Hamiltonian:  in  general,  the  charqes  at  each 
layer  do  not  coincide  with  the  initial  charges  shown  in  Figure 
2.  In  our  consistent  procedure,  we  allow  for  charges  in  the 
potentials,  VQ,  and  V 2»  Thus,  VQ,  and  V2  are  related  to 
the  charges  of  each  layer  through  a  kind  of  consistent  equation. 
This  means  that  is  defined  as  the  electrostatic  potential 
created  by  the  charges  Sn^,  as  measured  with  respect  to  the 
ideal  case  shown  in  Figure  2.  This  yields: 

V2  =  ad«n3 

=  ad  Sn^  +  ad(<5n2+6n3)  (1) 

Vq  =  ad6n3  +  ad(6n2+6n3)  +  ad 1  (  6n^+<5n2+  fin^  )  , 

where  d  is  the  interlayer  distance  in  GaAs  and  d'  is  the  distance 
between  the  Al  layer  and  the  last  As  layer;  a  =  4ir/A,  A  being 
the  area  per  surface  atom. 

Notice  that  in  our  actual  case 

6nQ  +  6n3  +  <$n2  +  6n3  =  0 

is  verified  with  a  very  good  approximation  since  the  interface 
perturbation  only  penetrates  three  layers  in  the  semiconductor. 

Eqs. (1)  define  our  consistent  conditions,  and  allows 
us  to  calculate  the  final  perturbations  V. . 


3.  RESULTS  AND  DISCUSSION 

Table  2  shows  the  values  of  VQ,  and  V2  as  obtained 

for  the  different  geometrical  configurations.  Figure  3  shows 

the  local  density  of  states  in  the  A1  layer  for  a  half  monolayer 

Figure  3a  shows  the  total  density  of  states,  while  Figures  3b 

and  3c  show  the  local  density  of  states  projected  onto  the 

p  +p  (bridge  bond)  and  the  s+p  (dangling  bond)  orbitals  of  Al. 
x  y  z 

From  these  figures  we  see  that  there  is  a  cation-like  surface 

band  in  the  semiconductor  gap;  this  band  is  a  dangling  bond 

state,  which  is  3/4th  occupied.  Notice  that  these  results  are 

similar  to  the  ones  obtained  for  a  cation-terminated  GaAs(lOO)- 
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surface  '  ;  the  change  of  Ga  by  Al  introduces  rather  small 

differences  between  the  two  cases;  in  particular,  we  find  the 
Fermi  level  pinned  by  a  high  density  of  states,  and  regions  of 
the  semiconductor  energy  gap  (outside  the  surface  band)  having 
a  zero  density  of  states. 

Figure  4  shows  the  local  density  of  states  in  the  Al 
atom  for  a  full  monolayer.  For  this  case,  and  comparing  with  a 
half  monolayer,  we  can  see  the  evolution  of  the  electronic 
density  of  states.  In  particular,  we  see  that  there  is  a  finite 
density  of  states  anywhere,  with  electronic  states  filling  the 
semiconductor  energy  gap;  however,  the  Fermi  level  for  this 
case  has  not  changed  very  much  from  the  half  monolayer  case: 
while  for  this  case  Ep  is  0.84  eV  above  the  valence  band  top, 
for  the  full  monolayer  we  find  that  Ef  is  at  1.15  eV.  On  the 
other  hand,  it  is  worth  stressing  that  results  found  in  other 
cases  23'2®  (ai  on  Si (111 )  and  Al  on  Ga(100))  show  that  the 


barrier  is  completely  formed  for  a  monolayer,  when  the  surface 


density  of  states  fills  completely  the  semiconductor  energy  gap. 
Based  on  these  results  and  in  our  present  calculations,  we  can 
conclude  that  the  Schottky  barrier  for  the  Al-GaAs(lOO)  junction 
is  completely  formed  for  the  full  monolayeri 

In  order  to  obtain  the  barrier  height  of  the  junction, 
we  note  that  our  calculations  yield  a  shift  of  the  Fermi  level, 
between  the  full  and  the  half  monolayer,  of  0.31  eV.  For  the  half 
monolayer,  our  results  also  show  that  the  Fermi  level  coincides 
with  the  one  obtained  for  a  cation-terminated  GaAs-(lOO)  surface. 
It  is  worth  stressing  that  our  calculations  yield  accurate  results 
for  the  shift  of  the  Fermi  level  when  depositing  the  metal  on  the 
semiconductor  14,15r  but  do  not  give  so  good  results  for  the 
initial  Fermi  level  (this  energy  is  very  much  dependent  on  the 
semiconductor  bulk  bands).  We  can  determine,  however,  the  Fermi 
level  at  the  interface  by  using  other  more  accurate  calculations 
for  the  free  cation-like  GaAs-(lOO)  surface.  In  Ref. 26,  it  is 
shown  that  the  Fermi  level  of  that  ideal  surface  is  around  0.55 
eV  above  the  valence  band  top.  Combining  this  result  with  the 
calculation  presented  in  this  paper,  we  conclude  that  the  Fermi 
level  at  the  Al-GaAs(lOO)  interface  is  located  at  0.86±0.05  eV 
above  the  valence  band  top,  yielding  a  barrier  height  of 
0 . 55i0 . 05  eV. 

Let  us  finally  comment  that  the  barrier  height  calculated 
in  this  paper  for  the  (100) -face  is  similar  to  the  one  calculated 
for  the  abrupt  Ag-GaAs(llO)  junction  30.  For  this  case,  the 
barrier  height  was  found  to  be  0.63±0.05  eV.  Since  changing  from 
Ag  to  Al  seems  to  yield  small  variations  in  the  barrier  height  14 r 
we  conclude  that,  to  an  accuracy  of  0.1  eV,  the  barrier  heights 


for  the  Al-GaAs(lOO)  and  Al-GaAs(llO)  are  the  same. 
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